Strong pressure-energy correlations in van der Waals liquids 
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Strong correlations between equilibrium fluctuations of the configurational parts of pressure and 
energy are found in computer simulations of the Lennard- Jones liquid and other simple liquids, but 
not for hydrogen-bonding liquids like methanol and water. The correlations, that are present also 
in the crystal and glass phases, reflect an effective inverse power-law repulsive potential dominating 
fluctuations, even at zero and slightly negative pressure. In experimental data for supercritical Ar- 
gon, the correlations are found to be approximately 96%. Consequences for viscous liquid dynamics 
are discussed. 
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For any macroscopic system thermal fluctuations are 
small and apparently insignificant. That the latter is not 
the case was pointed out by Einstein, who showed that 
for any system in equilibrium with its surroundings, the 
specific heat is determined by the magnitude of the en- 
ergy fluctuations. This result may be generalized, and it 
has long been well understood that linear-response quan- 
tities are determined by equilibrium fluctuations of suit- 
able quantities [l|, [2], yj • One expects few new insights to 
come from studies of fluctuations in equilibrated systems. 
We here report strong correlations between instantaneous 
pressure and energy equilibrium fluctuations in one of 
the most studied models in the history of computer sim- 
ulation, the Lennard- Jones liquid. These findings have 
significant consequences, in particular for the dynamics 
of highly viscous liquids. 

Using molecular dynamics 0, [H, fluctuations were 
studied for N — 864 particles interacting via the 
Lennard-Jones (LJ) pair potential [6j 4>Lj{r) — 
4e [(er/r) 12 - (cr/r) 6 ] in the NVT ensemble Q, where r is 
the distance between two particles. The configurational 
contribution to the instantaneous pressure defines the in- 
stantaneous virial W{t) by 0] p(t)V = Nk B T(t) + W(t). 
Fig. 1(a) shows normalized instantaneous equilibrium 
fluctuations of W(i) and the potential energy U(t) for 
a simulation at zero average pressure. The two quan- 
tities correlate strongly. To study the correlations sys- 
tematically, temperature was varied at five different den- 
sities. The results are summarized in Fig. 1(b), plot- 
ting instantaneous virial versus instantaneous potential 
energy, with each color representing equilibrium fluctu- 
ations at one particular temperature and density. The 
figure reveals strong W, U correlations with correlation 
coefficients mostly above 0.9, see table HI The results 
at a given density form approximate straight lines. The 
data include slightly negative pressure conditions, as well 
as three instances of the crystallized liquid (lower left cor- 
ner) . 

For any system with pair-wise interactions W(t) = 
-Y,i < j r ij(tW(r ij (t))/3%\$. Perfect W,U correlation 
applies if <fi(r) = ar~ n + <fio in which case AW(t) — 
(n/3)AU(t) where AW(t) = W(t) - (W), etc. An ob- 



vious first guess is therefore that the strong correlation 
directly reflects the r~ 12 term in the LJ potential. That 
is not correct because the exponent n = 12 implies a slope 
of 7 = 4 of the lines in Fig. 1(b); the observed slope is 
7 = 6 (±10%, see table HJ, corresponding to effective in- 
verse power-law exponents n w 18. The repulsive core of 
the LJ potential (r < 2 1 / 6 <r), however, can be well ap- 
proximated by <ppow (t) = ar~ n + c^o , with an exponent 
n considerably larger than 12 @, @|. If one requires that 
the O'th, l'st and 2'nd derivatives of the two potentials 
agree at r = r , one finds n(r ) = 6 + 12/[2 — (r /a) 6 ]. 
Thus, n(a) = 18, whereas n(0.969cr) = 16.2 (this is where 
4>lj = e). 

To directly test whether the fluctuations are well 
described by an inverse power-law potential, we pro- 
ceeded as follows. A large number of configurations 
from the simulation of the zero-pressure state-point in 
Fig. 1(a) were stored. This time-series of configura- 
tions was analyzed by splitting the potential energy into 
two terms: U(t) = U pow (t) + U rest (t), where U pow (t) = 
J2i<j 4'pow( r ij(t))i i- e - ; the potential energy if the inter- 
atomic potential were an inverse power-law. Comparing 
Up OW (t) and U(t) it was found that the fluctuations were 
nearly identical, AU pow (t) ~ AU(t), with a correlation 
coefficient of 0.94. Applying the same procedure to the 
virial, we found AW pow (t) » AW(t) with a correlation 
coefficient of 0.99. These results prove that the repulsive 
core of the LJ potential dominates fluctuations, even at 
zero and slightly negative pressure, and that at a given 
state point it is well-described by an inverse power-law 
potential. The fact that the repulsive forces dominate 
the physics - here the fluctuations - confirms the phi- 
losophy of the well-known Weeks, Chandler, Andersen 
approximation [Io| . 

It should be stressed that our approach is not to choose 
a particular inverse power-law and analyze the results in 
terms of it. In fact, for an exact inverse power-law poten- 
tial all data points in Fig. 1 (b) would fall on the same line 
[W(t) = ^U(t)]. Instead we simply study the equilibrium 
fluctuations at each state-point and find strong W,U cor- 
relations, which in turn can be explained by an effective 
inverse power-law dominating the fluctuations. The effec- 
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FIG. 1: Results from equilibrium molecular dynamics sim- 
ulations of 864 particles interacting via the Lennard-Jones 
potential studied in the NVT ensemble where "Argon-units" 
were used (a — 0.34nm, e — 0.997kJ/mol). (a) Normalized 
fluctuations at T = 80K and zero average pressure (den- 
sity = 34.6 mol/1) of the virial, W, AW(t)/ ^(Aff) 2 ) , and 
of the potential energy, AU(t)/y/ ((AU) 2 ). The equilibrium 
fluctuations of virial and potential energy are strongly cor- 
related, as quantified by the correlation coefficient: R = 
(AWAU)/^/((AW) 2 )((AU) 2 } = 0.94, where averages are 
over the full length of the simulation (10ns) after 10ns of equi- 
libration, (b) Configurational parts of pressure and energy - 
virial versus potential energy - for several state points of the 
Lennard-Jones liquid. Each color represents simulations at 
one particular state point where each data point marks in- 
stantaneous values of virial and potential energy from a 10ns 
simulation. The black dashed lines mark constant density 
paths with the highest density to the upper left (densities: 
39.8 mol/1, 37.4 mol/1, 36.0 mol/1, 34.6 mol/1, 32.6 mol/1). 
State points on the blue dashed line have zero average total 
pressure. The plot includes three crystallized samples (lower 
left corner). 



tive inverse power law exponent is weakly state-point de- 
pendent, and the above explanation is consistent with the 
qualitative trends seen in tableHJ Increasing temperature 
along an isochore or increasing density along an isotherm 
results in stronger correlation and smaller slopes, corre- 
sponding to a numerically smaller apparent exponent. 
This reflects particles approaching closer to each other, 
and thus rrj decreasing (n(rrj) decreasing) and the inverse 
power law being an even better approximation to the LJ 



p = 34.6 mol/1 


60K 


80K 


100K 


1000K 


R 


0.900 


0.939 


0.953 


0.997 


7 


6.53 


6.27 


6.08 


4.61 


T = 130K 


32.6 mol/1 


36.0 mol/1 


39.8 mol/1 




R 


0.945 


0.974 


0.987 




7 


6.06 


5.71 


5.40 




p = 0.0 GPa 


60K 


70K 


80K 


90K 


R 


0.965 


0.954 


0.939 


0.905 


7 


6.08 


6.17 


6.27 


6.52 



T ABLE I: The corre lation coefficient, R, and the slope 7 = 
\J ((AW) 2 } I ((AU) 2 } along an isochore, an isotherm, and an 
isobar for the Lennard-Jones Liquid. 



potential close to tq. We do find that 7 — » 4 (n — » 12) at 
high temperatures and/or densities as expected, but only 
under quite extreme conditions, see table |TJ Along an 
isobar there is competition between the effects of density 
and temperature. Our results show that the density ef- 
fect dominates: the correlation increases with decreasing 
temperature. This, incidentally, is the limit of interest 
when studying highly viscous liquids (see below). 

In order to investigate how general the W, U correla- 
tions are, several other systems were studied. If Wit) 
and U(t) are perfectly correlated (R = 1), the following 
identity applies: (AWAC/) 2 = ((AW) 2 ) ((AU) 2 ) . Fig. 2 
summarizes our simulations in a plot where the diagonal 
corresponds to perfect correlation and the y-variable by 
the fluctuation-dissipation theorem equals T times the 
configurational pressure coefficient [(T/V)(dW/dT) v \. 
Liquids with strong W, U correlations (R > 0.9) in- 
clude: 1) A liquid with exponential short-range repulsion; 

2) The Kob- Andersen binary Lennard-Jones liquid [ilT |: 

3) A liquid consisting of asymmetric "dumb-bell" type 
molecules (two unlike Lennard-Jones spheres connected 
by a rigid bond [I3|); 4) A 7-site united-atom model of 
toluene [l3| • The last three liquids are examples of good 
glass-formers that can be cooled to high viscosity without 
crystallizing. Liquids not showing strong W U correla- 
tions are methanol Q and SPC/E water [la ]: in these 
models the instantaneous potential energy has contribu- 
tions from both LJ interactions (ULj(t)) and Coulomb 
interactions (Uc{t)). Since the Coulomb interaction 
is an inverse power-law with n = 1, the correspond- 
ing contribution to the instantaneous virial is given by 
Wc(t) = Uc(t)/3, i.e., perfect correlation. For the LJ in- 
teraction of SPC/E water we find AW L j(t) « 6AU L j(t) 
with correlations coefficients above 0.9. Since the propor- 
tionality constants are different, however, the sums of the 
contributions do not correlate very well. In fact, close to 
the density maximum of water we find that W(t) and 
U(t) are uncorrelated. For methanol AW(t) and AU(t) 
correlate well at such high temperatures that the LJ in- 
teractions completely dominate («3000K). 
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FIG. 2: (AWAU)/(k B TV) plotted as a function of 
({(AW) 2 )((AU) 2 )) 1/2 /(k B TV) for several liquids. If the cor- 
relation is perfect (R = 1) the data fall on the diagonal. 
The region above the diagonal corresponds to R > 1 and 
is thus forbidden. LJ: Lennard- Jones results from the sim- 
ulations reported in Fig. 1. Similar results were found in 
the NVE ensemble, for larger samples, and using Langevin 
dynamics (results not shown). Exp: 500 particles interact- 
ing via a pair potential with exponential repulsion; U(r) — 

| |6e~ 14(r/<T " 1) - 14(cr/V) 6 ] , simulated with Metropolis dy- 
namics. Kob- Andersen binary LJ: The Kob- Andersen 
binary Lennard- Jones liquid, N=1000 [llj. This includes 
data for the less-viscous liquid, the highly viscous liquid, as 
well as the glass. Asymmetric dumb-bell: 512 asymmet- 
ric "dumb-bell" molecules Q3]. OPLS-UA Toluene: A 7- 
site united-atom model of toluene, N=1000 [II]. GROMOS 
Methanol: 512 methanol molecules Q. SPC/E water: 
4142 SPC/E water molecules [l]|. Except for the two systems 
with Coulomb interactions (Methanol and water), all systems 
studied have strong correlations between fluctuations in virial 
and potential energy; correlation coefficients are above 0.9 for 
all state-points shown, except those with negative pressure, 
where they are slightly smaller. The correlation coefficients 
increase with increasing density and temperature. 



Do strong pressure-energy correlations have conse- 
quences accessible by experiment? In the following we 
demonstrate how it is possible to test for strong W,U 
correlations in systems where the kinetic contribution to 
the isochoric heat capacity is known, exemplified by ex- 
perimental data for supercritical argon. From the defi- 
nition of R, utilizing three fluctuation formulas Q it is 



straightforward to show [16| (where c' v 



,conf 



and (3 



conf 



the configurational parts of the isochoric heat capacity 
per volume and pressure coefficient respectively, and Kt 
is the isothermal bulk modulus) that: 



FIG. 3: Data for supercritical Argon at 3 different densi- 
ties covering the temperature range 200K-660K jUJ show- 
ing &W,U correlation of 96% [K T = -V{dp/dV) T , p coni = 
p - Nk B T/V = W/V, /^ onf = (dp/dT)v - Nk B /V, and 
c^ ni = C v /V - (3/2)Nk B /V\. The full line corresponds to 
perfect correlation between virial and potential energy [lr3 ]. 
The inset shows the prediction for an exact inverse power-law 
interatomic potential (full line): B(T) = A(T). The poor fit 
shows that the fact that fluctuations are well described by an 
effective inverse power-law does not imply that this is the case 
also for the equation of state. 



T(^ onf ) 2 /c c v onf = R 2 (P-K T + (X)/V) 



(1) 



Here X = J2i<j rw'(ry)/9 is the so-called "hypervirial" 
where w(r) — r<p'(r). This quantity cannot be deter- 
mined experimentally Q , so we apply an approximation. 
For an exact power-law potential, one has X — (n/3)W, 
which, however, is not expected to be a good approxima- 
tion since the apparent power-law depends on the state- 
point. In the vicinity of a given reference state-point, 
however, one expects [X—X re f] ss (n/3)[W— W re f] along 
an isochore (confirmed by our simulations). Using this 
approximation and assuming that R is roughly constant, 
we can test for strong W,U correlations. Fig. 3 shows 
experimental data for supercritical Argon covering the 
temperature range 200K-660K at three different densities 
17| . showing that W and U correlate 96% in this case. 



The apparent power-law exponents, n = 3/3y nf /c™ nf jl6| 
varies from 13.2 to 15.8, decreasing with increasing tem- 
perature and density. The inset shows that the Argon 
data do not follow the prediction following from an exact 
inverse power- law potential [3, [T3| (Eq.© with R = 1 
and X = (n/3)W); thus the W,U correlations show that 
such an effective power-law description applies to a good 
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approximation only for the fluctuations. This is analo- 
gous to the situation for the Lennard- Jones liquid: The 
equation of state is poorly described by that following 
from an inverse power-law potential (see e.g. [Ill), al- 
though the fluctuations are well described by this. 

A different class of systems where it is possible to test 
for strong WJJ correlations experimentally is highly vis- 
cous liquids H, 20| . These are characterized by a clear 
separation of time scales between the fast vibrational de- 
grees of freedom on the picosecond time scale and the 
much slower configurational degrees of freedom on the 
second or hour time scale, depending on temperature 
Hi HI 0, HI HI US W\ . Suppose a highly viscous 



liquid has perfectly correlated W, U fluctuations. When 
W and U are time-averaged over, say, one tenth of the liq- 
uid relaxation time [13], they still correlate 100%. Since 
the kinetic contribution to pressure is fast, the time- 
averaged pressure equals the time- average of W/V plus a 
constant. Similarly, the time- averaged energy equals the 
time-averaged potential energy plus a constant. Thus 
the fluctuations of the time-averaged W and U are the 
slowly fluctuating parts of pressure and energy, so these 
slow parts will also correlate 100% in their fluctuations. 
This is the single "order" parameter scenario of Ref. 
[20j ]. In this case, knowledge of just one of the eight 
fundamental frequency-dependent thermoviscoelastic re- 
sponse functions implies knowledge of them all [2(1 (ex- 
cept for additive constants [Hi)- This constitutes a con- 
siderably simplification of the physics of glass-forming 
liquids. Unfortunately, there are few reliable data for the 
frequency-dependent thermoviscoelastic response func- 
tions 29]. Based on the results presented above we 
predict the existence of a class of "strongly correlating 
liquids" where just one frequency-dependent thermovis- 
coelastic response function basically determines all. Our 
simulations suggest that the class of strongly correlating 
liquids includes van der Waals liquids, but not network 
liquids like water or silica. This is consistent with the 
findings of De Michele et. al. [30j |. 

Very recently Coslovich and Roland studied diffu- 
sion constants D in highly viscous binary Lennard- 
Jones mixtures at varying pressure and temperature 
3l| . Their data follow the "density scaling" expression 
D = F{p 1 /T) [32j, and they showed convincingly that 
the exponent 7 reflects the effective inverse power law of 
the repulsive core. In view of these findings, we conjec- 
ture that strongly correlating viscous liquids obey density 
scaling, and vice versa. If this conjecture is confirmed, 
by virtue of their simplicity the class of strongly correlat- 
ing liquids provides an obvious starting point for future 
theoretical works on the highly viscous liquid state. 
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